Dimensional trend in CePt2ln7, Ce-115 compounds, and Celns 
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We present realistic Kondo-lattice simulation results for the recently-discovered heavy-fermion 
antiferromagnet CePt2ln7 comparing with its three-dimensional counterpart Celns and the less 
two-dimensional ones, Ce-115's. We find that the distance to the magnetic quantum critical point is 
the largest for Celns and the smallest for Ce-115's, and CePt2ln7 falls in between. We argue that the 
trend in quasi-two-dimensional materials stems from the frequency dependence of the hybridization 
between Cerium 4/-electrons and the conduction bands. 

PACS numbers: 71.27.-|-a, 74.10.-l-v, 75.40.Mg 



Recently-discovered CePt2ln7 [l|, Q has provided a 
new approach toward the two-dimensional (2D) limit 
in Celna-derived heavy-fermion material family, among 
which Ce-115 materials have been discussed intensively 
in the past decade. In the hot debates on 115's, the 
possible scaling of the superconducting transition tem- 
perature Tc to c/a, where a and c are the lattice con- 
stants of the tetragonal crystal structure, has been dis- 
cussed [1,13] and there has been a hope to raise Tc by mak- 
ing a more 2D-like material. In this context c/a = 4.694 
(if c should be taken to be the interlayer distance be- 
tween Cerium planes, the scaling parameter should be 
taken to be 2.347) in CePt2ln7 |3j looks promising as 
compared to the typical values of c/a ^ 1.6 in Ce-115's 
where the highest T^s have been found among Cerium 
heavy-fermion compounds [^-01 ■ 

Theoretically, the superconductivity (SC) with high 
Tc's in strongly-correlated materials has been discussed 
to be mediated by magnetic fluctuations [8|] and the pos- 
sible mechanism for SC seems to be intimately related to 
the nearby magnetic quantum critical point (QCP) 
An empirical material-designing principle to get high Tc 
would be to make it as close as to QCP and it is desir- 
able for a given possibly high-Tc material to know how 
close it is located to QCP. The robustness of the mag- 
netic pairing has been discussed considering the spatial 
dimensional effects IQ], stressing the stronger robustness 
in the more 2D-like systems even though care must be 
taken regarding the details of the electronic structure. 
Thus we are motivated to address QCP in CePt2ln7 fully 
taking into account its electronic structure and see the 
trend among the related materials, namely, Celns with 
c/a = 1 and Ce-115's with slightly larger c/a ~ 1.6, to 
see if any microscopic origin of c/a-scaling of Tc could 
be traced quantitatively to QCP. In the present work we 
predict the QCP's for all of these materials to elucidate 
where exactly CePt2ln7 is located in the neighborhood 
of QCP. 



Recent experiments on CePt2ln7 [llj have shown that 
this material is an ideal Kondo lattice with commensu- 
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FIG. 1: (Color online) Hybridization function calculated by 
local density approximation (LDA) with Hubbard-1 approxi- 
mation for CePt2ln7, Ce-115's and Celus using the code de- 
scribed in Ref. [31 . The trace in the vertical axis quantity is 
taken over 14 orbitals of the 4/-orbital of Cerium. 



rate antiferromagnetism. It can be a good target for 
realistic Kondo lattice simulations [13] for which it has 
been shown that the description works fine in the Kondo 
limit, meaning that local /-level position, ef, measured 
from the Fermi level should be negative and large, U + ef 
should be positive and large with U being the on-site 
Coulomb repulsion energy, and at the same time hy- 
bridization should not be too big. When the va- 
lence fluctuations start to dominate the Kondo lattice 
simulations does not work and we have to go back to 
the lattice simulations based on the Anderson impurity 
problem [11], which would increase the computational 
cost. It has been discussed that Ce-115's, especially 
CeRhIn5 and Celrlns, have well localized electrons |14l ] 
and we expect that the comparison between CePt2ln7, 
CeRhlns/Celrlns, and Celns in our simulations would 
make sense, while CeCoIns might have to be looked at 
with some extra care. 

We define our realistic Kondo lattice Hamiltonian 
starting with the first-principles electronic-structure cal- 
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culation for the Cerium heavy fermion materials. The 
local density approximation (LDA) for delocalized s-, p-, 
and d-conduction electrons combined with the Hubbard- 
I approximation ^J3\ for localized /-electrons gives the 
frequency-dependent hybridization function, — QA(a;)/7r, 
between the conduction electrons and the localized 4/- 
electrons as shown in Fig. [1] A part of the trend is al- 
ready seen at this stage: we find that the hybridization 
around the Fermi level, which is most relevant for Kondo 
screening is the strongest for Celns and the weak- 
est for CePt2ln7, and Ce-115's comes in between. This 
is in line with the relevance of the off-plane hybridiza- 
tion that has been discussed for Celrln5 [3]. Thus ob- 
tained data reflects the high-energy physics of the given 
material and then we obtain our low-energy Hamilto- 
nian via the Schriefler- Wolff transformation imple- 
mented in a realistic way fl^]. The Kondo coupling 
Jk = y^[l/|e/| + l/{Uctf + £/)] incorporates the vir- 
tual hopping process — > /° in the first term and in 
the second term another virtual process with 
Uctf = U — JhuiicI ^ 4 [eV] being the effective on-site 
Coulomb energy in the multiplet in the Ce atom in- 
cluding the effective Hund coupling Jnund ^ 1 [eV]. The 
hybridization is taken to be/_^^ [— 9A(w)/7r], to 
take care of the normalization of the density of states 
of our conduction electrons with keeping the value of a 
relevant quantity in Kondo physics, Jkp(w), where p{uj) 
is the conduction electron density of states, to the given 
realistic value. Note that the conduction band cutoff D 
matters in the definition of Jk- Although — 3A(a;)/7r is 
a decaying function of 0, it does not fall off suffi- 
ciently fast due to numerically imperfect projection onto 
/-states. Thus the cutoff D is implemented to be equal 
to the on-site Coulomb energy D = U = 5 [eV] . The rele- 
vant parameters and the outputs of LDA for the Cerium 
compounds of our present interest is summarized in Ta- 
ble m We note that our calculated values of for Ce- 
115's are not exactly consistent with the photoemission 
data in the literature 1 



26j . However both of the calcu- 
lated values and experimental ones indicate that all of the 
materials are deep in the Kondo regime. The spin-orbit 
splittings are set to be 0.3 [eV] for all of the materials. 
The structural parameters to be plugged into LDA are 
taken from experiments referring to Ref. (25l |. 

Now we describe how we solve the realistic Kondo lat- 
tice model (KLM) that we have defined from — 3A(a;)/7r 
in Fig. [1] and Table HI We use dynamical mean- field the- 
ory (DMFT) [27, Q that is formulated on a local /- 
electron basis [29| which enables us to reach low temper- 
ature region with a modest computational cost and ad- 
dress the QCP in a semi-quantitative way utilizing state- 
of-the-art continuous-time quantum Monte Carlo (CT- 
QMC) impurity solver sO-ll]- We plug- in the realistic 



TABLE L The parameters for the investigated materials. The 
first two columns are LDA -I- Hubbard-I results for the local 
/-level position e / and the value of the hybridization function 
on the Fermi level. The final column shows the crystal-field 
splittings in the tetragonal structure. 
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CePtalnr 


-1.81 


0.174 


8.6, 12.9° 




CeCoIns 


-1.97 


0.205 


6.8, 25'' 




CeRhlns 


-1.90 


0.209 


5.9, 28.5'= 




Celrlns 


-1.95 


0.220 


5.26, 25.875^= 


3D 


Celna 


-1.72 


0.239 


12"^ 



" = lOOfKJand 150 [K], a rough estimate 

'' Ref. [2l|,[2i|. Another crystal-field scheme in Ref. ^ gives 

essentially the same results. 

Ref. [24| . The latest crystal- field schemes in Ref. ]21j give 
the values close to these. 

Ref. [2^. This materials has the cubic structure which 
brings Ai — A2. 



of DMFT. We look at the magnetic phase transition as 
a function of temperature and restore the Doniach phase 
diagram [34*1 by varying the Kondo coupling Jk for a 
given material. Thus a magnetic QCP is found on a 
realistic Doniach phase diagram spanned by the Kondo 
coupling and temperature, and the realistic data point 
for the given material is picked up for the realistic value 
of U — 5 [eV] and JhuikI = 1 [eV] to estimate its distance 
to QCP. 

Now we show how we determine the QCP of CePt2ln7. 
The calculated temperature dependence of the inverse 
of the staggered magnetic susceptibility, l/x(7r), for 
CePt2ln7 is shown in Fig. [2l Here we have employed 
random dispersion approximation 35| to estimate the 



crystal-field and spin-orbit splittings as given in Table H] 
in the local 4/-level in the impurity problem. Thus our 
solutions are numerically exact up to the approximation 



two-particle Green's function by decoupling, which would 
enhance our transition temperatures on top of DMFT as 
will be seen below. Being consistent with Doniach's pic- 
ture (S^l where the winner of the competition between the 
magnetic-ordering energy scale oc J^p and the Kondo- 
screening energy scale oc exp[— 1/( Jkp)] interchanges at 
some finite Jk, it is seen that small Jk's give a diverg- 
ing x(7r) at a finite Neel temperature Tn while large Jk's 
give a saturating xC"") low temperatures, and some 
value of Jk in between gives the quantum critical point 
where Tn vanishes. The Neel temperatures for smaller 
Jk's are determined by linear extrapolation of l/x(7r) to 
the lowest temperature region and thus obtained Tn is 
plotted against the Kondo couphng in Fig. [3] in a format 
of restored Doniach phase diagram. Here we note that 
we just identified what we call QCP by a parameter seg- 
ment where the finite Neel temperature seems to have 
vanished, and there is always a possibility that in some 
smaller parameter segments there is actually a coexis- 
tence region. There is also a possibility numerically for a 
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FIG. 2: (Color online) Temperature dependence of the in- 
verse of the staggered susceptibility in the realistic Kondo lat- 
tice models for CePt2ln7 by which we identify that the phase 
boundary is located in 1.33Jk,o < Jk < 1.34Jk,o, where Jk,o 
is the Kondo coupling which corresponds to U — 5 [eV] and 
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FIG. 3: (Color online) Restored realistic Doniach phase dia- 
gram for CePt2ln7 from the data in Fig. (2] AF is the anti- 
ferromagnetic phase and PM is the paramagnetic phase. The 
line is a guide to the eye. 



first-order phase transition where the Neel temperature 
actually jumps from a finite value to zero at a certain 
point of Jk- We leave the exact characterization of what 
we also call QCP here for future investigations and for 
the moment we would be satisfied with that it looks like 
QCP practically in most cases with a very small jump 
numerically, if any. 

For CePt2ln7, the realistic data point is at Jk = 
1.165Jk,0) which corresponds to U — 5 [eV] and Jnund = 
1 [eV], gives the Neel temperature Tn — 10 [K] which is 
larger than the experimental result Tn — 4.5 [K] 11] due 
to our mean-field argument with DMFT and the random 
dispersion approximation 35[ which becomes exact on a 
lattice with perfect nesting. The same analyses are ap- 
plied to the other materials to give the Neel temperature 
Tn = 19 [K] for Celng, 15 [K] for CeRhlns, and 14 [K] for 
CeCoIns. Again the experimental values (Tn = 10.2 [K] 
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FIG. 4: (Color online) Realistic Doniach phase diagram for 
CePt2ln7, and Ce-115's together with their parent compound, 
Celns with the horizontal axis rescaled to measure the dis- 
tance to the magnetic QCP. The line is a guide to the eye. 



for Celna, 3.8 [K] for CeRhlng 

comes below these 

results. The data for CeCoIns is not consistent with the 
experimental fact that CeCoIns is a non-magnetic heavy 
fermion material @. However the arrangement of the 
materials around QCP reveals that the apparent finite 
Tn actually comes from being in an immediate proximity 
to QCP as shown in Fig.Sl We use a rescaled Kondo cou- 
pling on the horizontal axis in terms of its value right on 
the QCP, t = (Jk - Jk|qcp)/>/k|qcp, to remove a prob- 
lem with the Doniach phase diagram that each material 
has its own energy scales. 

Now we inspect the distance to the QCP of CePt2ln7 
referring to those of Ce-115's and Celna. The cubic par- 
ent material Celn^ is seen to be most separated from 
QCP and Ce-115's are found to be concentrated in the 
neighborhood to QCP, with CeRhlns on the magnetic 
side and Celrlns on the non-magnetic side. Our numeri- 
cal resolution is not sufficient to locate CeCoIns in its cor- 
rect non-magnetic side, but it is clear that it sufficiently 
works to estimate the extreme closeness to QCP of Ce- 
115's. We note that our calculation scheme might not be 
as good for CeCoIns as for the others due to the possibly 
stronger effects of valence fluctuations in this material. 
It is seen that when we try to reach the quantum critical 
point (QCP) from Celn^ in the three-dimensional limit, 
CePt2ln7 is located in the midway toward Ce-115's which 
are located closest to the QCP. Seen from QCP, Celna 
is already close enough to enable the pressure-driven su- 
perconductivity [36| and CePt2ln7 would also have one. 
Making a material more 2D indeed helps it to come closer 
to QCP, which is reasonable in the general context that 
the lower spatial dimensionality would suppress the mag- 
netic long-range order. However within the 2D-side, the 
trend among CePt2ln7 and Ce-115's is somewhat non- 
monotonic. The possible reason is seen in Fig. [1] where 
a dip around the Fermi level is seen for CePt2ln7 which 
would help to drive it to the magnetic side with the re- 
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duced Kondo screening. On top of the spatial dimen- 
sionality, the frequency dependence of the hybridization 
seems to introduce the nontrivial trend in this way. 

Here we note that precisely speaking we are discussing 
in the infinite-dimensional limit with the DMFT, but at 
least a semi-quantitative trend of Tn would be satisfacto- 
rily addressed. We also note that in the literature the im- 
portance of momentum dependence of the hybridization 
has been stressed 3 37, 38| which we have neglected. 
What we have seen is that the energy dependence seems 
to be sufficient at least to capture the trend in the dis- 
tance to QCP, thus for the prediction of magnetically- 
mediated superconductivity with a possible high Te- 
la order to reach QCP, it would be interesting to have 
a material analogous to CePt2ln7 but without a big dip 
in the hybridization around the Fermi level. For that, the 
following possible ways for the material designing could 
help: 1) electronic carrier doping, 2) ascending T in the 
periodic table for CeT2ln7 to enhance the hybridization, 
and 3) shifting T to the left- or the right-hand side on 
the periodic table to lift the dip off the Fermi level to 
enable the stronger Kondo coupling to drive the material 
toward the QCP side. 

Finally we discuss the degree of the delocalization of 
/-electrons by looking at their contribution to the Fermi 
surface of the conduction electrons in our realistic Kondo 
lattice description for each material. Even if we have only 
localized /-electrons in our Hamiltonian, they contribute 
to the formation of the large Fermi surface via the hy- 
bridization and in this sense they can be delocalized [39j . 
We follow the procedure used in Ref. [13] ■ From our sim- 
ulations for a give material at a fixed temperature, we 
get the conduction-electron self energy E(ia;„) and then 
track the temperature dependence of the real part of it 
at iuJn = to look at the shift of the Fermi level. Here 
ujn = (2n + 1)t:T is the fermion Matsubara frequency. 
The result is plotted in Fig. [S] It is seen that CePt2ln7 
has the most localized /-electrons down to the tempera- 
ture range of 20 [K] among these materials while the 115's 
have the strongest delocalization. Celua falls in between. 
We note that the trend in the /-electron delocalization 
does not exactly follow that in the hybridization strength 
right on the Fermi level seen in Fig. [T] and it is the out- 
come of the frequency dependence of the hybridization 
that we have taken into account in our realistic Kondo 
lattice simulations. 

To conclude, we have predicted the QCP of CePt2ln7 
and discussed the possible strategy to have the higher Tc 
in the related material family within our realistic Kondo 
lattice description which has been shown to predict the 
properties of Ce-115's and their parent material Celua 
semi-quantitatively. Our method's predictability of mag- 
netic QCP would further provide the material designing 
principle toward more high-Tc materials in the upcoming 
material exploration. 
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FIG. 5: (Color online) Temperature dependence of the real 
part of the conduction-electron self-energy for CePt2ln7, 
Cerium-115's and Celna. 
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